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ABSTRACT 


A differential equation for the scale tensor in turbulent flow 
is developed from basic considerations and applied to the flow of 
a constant-density fluid in the boundary layer on a flat plate. 
Results from preliminary runs of a computer implementation are 
discussed. 
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1.- INTRODUCTION 


The method of analysis of turbulent flows known as full second-order 
closure modeling has been applied with great success to a variety, 
of flows by A. R. A. P. and others. However, one weak point has 
remained in this method; that is, the determination of the scale 
which is an important factor in many of the models. Much, of the 
work has used an empirical function to determine the scale. 

Although this has worked well in many cases, it requires reevalua- 
tion whenever a novel situation is considered. 

Another approach has been to use a differential equation to follow 
the evolution of the scale. This too has been reasonably successful 
in many cases, but it requires considerably more modeling (and 
evaluation of modeling parameters) and the results achieved have not 
been such that universal validity is indicated. 

One way of arriving at such a differential equation starts by 
considering a two-point velocity correlation. This leads naturally 
to the consideration of a scale tensor, measuring the correlation 
length for the various pairs of fluctuations of velocity components. 
It is the contraction (or trace) of this tensor that serves as a 
measure of the scale in the modeling. 

It is evident that much of the information contained in the tensor 
is lost in the contraction process. During the past year, a 
computer program to solve the differential equations for the whole 
scale tensor has been developed at A.R.A.P. for Ames Research Center 
of NASA under Contract NAS2-8014, as an extension of an effort to 
develop and understand turbulence modeling, especially as applied to 
compressible boundary layers. It is hoped that an analysis of the 
results of this program will lead to a better understanding of the 
concept of scale in turbulence modeling and of the turbulence 
itself. In order to keep this first effort in the area at a 
reasonable level,, the project has been restricted to the flow of a 
constant -density fluid in the boundary layer on a flat plate. This, 
of course, has the advantage that such flows are well understood 
experimentally. 

Much of the theory upon which the program is based was developed 
with the support of the Air Force Office of Scientific Research 
under Contract F44620-76-C-0048. 

This paper presents a derivation of the equations that are integrated 
by the program, with a description of the modeling entailed, and 
includes a description of the results of some preliminary runs. A 
companion paper, reference 1, serves as a program manual with 
descriptions of the numerical methods, operating instructions, and 
so forth. In an attempt to make each paper self-contained, a 
certain amount of material appears in both. 
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General tensor notation is used in much of the analysis. This is 
done to emphasise that the modeling is intended to apply in any 
(constant-density) flow situation using any coordinate system. It 
is valid to think of the equations as Cartesian, with no distinc- 
tion between subscripts and superscripts , with all derivatives 
interpreted as simple partial derivatives, and with the metric 
tensor, , interpreted as the Kronecker delta, 6^ . (Note 

that 6^ * 3 , by the summation convention.) The subscript t 

denotes a derivative with respect to time. 


The next section describes the modeling used in the equation for 
the Reynolds stress, and the following section treats the derivation 
of the scale equation. In section 4, the reduction of these equa- 
tions to the case of the boundary layer on a flat plate is 
presented. Finally, section 5 discusses the evaluation of the 
modeling parameters and the results of the computer calculations 
made so far. 



It is a pleasure to acknowledge the help of Morris Rjubesin of NASA 
Ames Research Center. His encouragement and support over the course 
of several years have been instrumental in carrying out this 
project. Me also thank our colleague, Ashok Varma, who made signi- 
ficant contributions to several aspects of the analysis. 




2.- THE BASIC MODELS 

For a fluid of constant density, the equations of notion are 
Continuity 


- ° 


(l) 


Momentum 


£ Pi 

\ + U U i,£ = 


+ VU.’* 

p ^ 


( 2 ) 


In this instance, it is convenient to rewrite equation (2) as 


u. + (u u.) = - P . + vu. 




1 »A 


, 3 - 


i 


(3) 


where equation (1) has been invoked to modify the second term and 
P has been substituted for p/p . P , sometimes called the 
"kinematic pressure," is loosely called just the "pressure" in what 
follows. The variables are written as sums of means f denoted by 
bars, and fluctuations, denoted by primes: 


u. * u. + u! 

i i i 


P = P + P* 


(4) 

(5) 


The fluctuations then are the deviations from the means, and by t 
definition the mean of a fluctuation is zero. It is assumed that 
the operation of taking the mean or averaging is interchangeable 
with differentiation. 

The substitution of equations (4) and (5) into equations (1) and 
(3) gives 


-£ , , £ « 
U ,£ + U , £ = 0 


( 6 ) 


r/)t 

iu* 


iu 


tr.tr. 


\ * x + + w'H'.t + 




vt 

t * «« 1 


p ,i - p 'i + vi i'!» + vu i’!i 


( 7 ) 


Taking the mean of these two equations results in 


*\i ~ 0 


( 8 ) 


X + + (u ’ Tu i ) ,£ 


- p ,i + v v!* 


(9) 


Equation (8) is identical in form to equation (1) , but equation (9) 
has a term with no counterpart in equation (3) . This term, involv- 
ing the correlation of fluctuations, is often called the turbulent 
stress or Reynolds stress. 

For engineering purposes, the main interest lies in determining the 
mean flow; to do so involves determining the Reynolds stress. 

Until about ten years ago, the only practical approach was to model 
the extra terra in equation (9) in terms of the mean flow (using 
Prandtl's "mixing length" concept, or "eddy viscosity")- This 
technique, which is still very useful in many circumstances, is 
cdlled first-order closure. However, it is recognized that in some 
circumstances this is not, indeed, cannot be, sufficient. The next 
step then is to derive the equation that determines the evolution 
of the Reynolds stress. 

o' * V <• 0 .' p o - I : .. o o t> ' c . •. > . o; C , r r 'v ; ..vo- r. ° C •; . ’ . 

Subtracting equations (8) and (9) from equations (6) and (7) gives 

«’% = 0 ( 10 ) 


u‘ t + Cu*up f£ + 







+ vu! 




(11) 


Rewrite equation (11) with j substituted for i and multiply by 
u| . Add the result to the original equation (11) multiplied by 

uj . Making use of equations (8) and (10) and averaging, the 

result is 


**£ / i t 


. * * 


+ + U 'X G j,t + u ’*“j 5 i,t + C»'X u p|i 




- u i p :, - “i p :i + • . ■» + ) <12> 


Using the rule for the derivative of a product, the right-hand side 
of this equation can be rearranged to give 


(u!u!) t + ^(uluj)^ + u’*u’ u j>£ + u**uj £ + (u^ujuj)^ 


,£ , - 


* * 


i 9* % t < 


(P‘uj) . 

x I J 


(p ‘ u i ) .i + p,u M + p,u i.i 


+ v(u!u-)-V 2v u! c u!’ £ 


ij 


j 


(13) 


An attempt to find an equation for 


U i U j 


has introduced several 


more unknowns , a triple correlation, correlations involving P’ 
and a correlation of derivatives of the velocity fluctuations . 
Second-order closure consists of devising models for these terms 
in such a way that the new equation and equations (8) and (9) form 
a closed set. v 'a .„<* 


In order to ensure applicability over a wide range of circumstances, 
the models used are invariant under Galilean transformations and 
vender coordinate transformations, they respect symmetry and, of 
course, they are dimensionally consistent. The "realizability" 
conditions recently pointed out by Schumann (ref. 2) have not been 
checked for these models except partially and incidentally. It 
appears from Schumann’s results that a detailed check might provide 
some useful constraints on the values of the parameters. ■ 4 . . 

Two quantities, A and q , appear in most of the models. The 
determination of the first of these, A , the scale, is the motiva- 
tion for this whole study. The square of the other, q , is 
defined by 


|*<V 
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so it is a measure of the energy of the turbulence. 

Most of the modeling constants are designated by rather cumbersome 
three -or four-letter names suitable for direct use in FORTRAN. 

The particular names used here were arrived at through an attempt 
to systematize the naming of the much larger set of constants 
needed in the compressible case (ref. 3). 

The model for triple correlations is practically the same as that 
given in reference 4: 


’ - VD0 <i A (< u i u j> ,k + j. + W> ,j) <14) 

Note that the three terms are needed for symmetry. 

The model for the dissipation term given in reference (4) is 

r u!u! 

u! ,u! ’ «* — (15) 

1 ’ x J 

where X is another scale given by 

,2 A 2 

= a + b M; 

te ( a 0 and b being constants) Thus , equation ' (15) can ..be ' written 


u! 


u* ’ ^ =* 


a 

tr 


ST5T 




. u -t u i 

Av i j 


(16) 


As was mentioned in reference 4, an isotr opic dissipation model in 
which g^.q 2 /3 is substituted for uTuT in equations (15) and 

(16) , had been tried and found to lead to negative values of autocor 
relations in some circumstances. (That is, we had violated one of 
the realizability conditions.) Eventually it was realized that the 
arguments favoring isotropic dissipation apply only at high values 
of 


6 


( 17 ) 


and that the problems had occurred only at low values, so g^j<l 2 /3 

was substituted for u.ju! in the second term on the right in 

equation (16) but not in the first. Actually, the formula in use 
allows a choice through another parameter, 6 : 


u! -u! ’ A ** -y u!u! + 

J a 2 1 3 Av 


“i“j + 


1 - 6 


; lj q2 ) 


Thus, 6=1 corresponds to equation (15) or (16) and 6-0 to 
isotropic dissipation in the high Reynolds number term only. 

The models involving the pressure fluctuation have been expanded 
considerably from those given in reference 4. They have been 
suggested by the following considerations. If the divergence of 
equation (11) is taken, many terms drop out by virtue of equation 
(10). Solving the resulting expression for the Laplacian of P 
and using equation (8) gives 


£ * 


- ,n - u ,£ u ,n 0 + (u ,£ u ,n ) 


- 2u u 


,n u ,£ ■ ,n ,£ 


Regarding this as Poisson’s equation for P' , the solution can be 
represented, in the absence of wall effects, as an integral 


P’(x) « - 



x - £| f ^ 




where S' represents the right-hand side of equation (19). The 
correlation of P' with another fluctuation, say <J>' , can thus be 

written : ..vr'..- •. 


* ,p ‘ - ' z^r 


il(i) s^(p d | 
|x - cl 


Recognizing that the two-point correlation under the integral can 
be expected to be negligible outside a region of characteristic 
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size A , and that the value of S' at £ * x can be assumed to 
be representative, this expression is written 


*'P’ - CA 2 0'S’ « CA 2 ( 2u £ $'u ,n 0 + $*u ,A u ,n e ) 

where C_ is a dimensionless quantity of order one. The P’ 
correlations needed for equation (13) are, therefore, written 

FF = A 2 (PMU u £ n u!u ,n 0 + PMU2 u!u’* u ,n „ ) 

■** v x x ,n / 


( 20 ) 


P ' U i.j * a2 ( PGD “\n + PGU2 < 21 > 


The four constants introduced here should, by the arguments just 
given, be all equal to C or 2C . They have been given separate 
names in the interest of flexibility and in particular to allow 
studies of the effect of separate terms . 

It appears that not much progress lias been made; in modeling two 
correlations, four new ones that need modeling have been introduced 
However, the nature of some of the new ones is such that the model- 
ing is quite constrained. This is particularly true of uTuT~T . 
When i and j are equal 1 


u i u i,k “ 2 (u i u p,k (not SUTnmed) 

"must: hold^ On the other hand, if j and k are contracted, 


.. 


u i u ’ % * 0 

by equation (10). Putting these conditions together practically 
forces 


“L u j,k - I K“j>, k - s hk<“ ,<! “i>;i - (22) 


a model remarkable for having no adjustable constants. 
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Similar considerations lead to 


u i,j u n,m ( ®mj u i u n * 3 g nm U i U j ’ 3 8 ij u n u m + S' q ) 


WGUB 


f- (gmj u i u ; * V U j U m) 


WGUC / „ n 2 _ 1 

72 \®rai®nj q 3 ®ij 


s«s 2 ) 


(23) 


The next has just one free index, so it is modeled by 


- « HU1 !( u ’ lu i).it 

Finally, the model 


(24) 


u! .u ,£ u ,n 0 = WWGU 
i.J ,n ,1 


M“i”i q2 ) <25) 


is governed by the requirement that it be zero when i and j are 
contracted. 

Substitution of equations (22) through (25) into (20) and (21) with 
the use of equation (8) gives 


p ‘ u i -W? a2 (| ^,n (u i u * n) ^ + \ 


'>’•0 £** C &~Ji- . o ° o>» 

+ PUW 


.. « » .. . £9> %' «• 


q A ( u * £ u!) t , 


v <: V v t ’ V 

(26) 


P’u! . = GUA 
1 » J 


(“j.n “i U ' n ‘ 3 


+ SUB (ij „ u!u' n - i* (i »]»i)+ SUC S M , 2 
+ I SU2 f(uluj - j g<y q 2 ) <27) 




/" 


where PUW - PMU2 .. WWU1 

GUA = PGU •. WGU 
GUB = PGU • WGUB 
GUC * PGU • WGUC 
GU2 =2 * PGU2 • WWGU 


Equations (26) and (27) are taken to be models. Hie steps leading 
to them are considered to be merely suggestive. The last term in 
each is essentially the corresponding model in reference 4. The 
GUA, GUB, and GUC terms in equation (27) are equivalent to terms 
arrived at by a different route by Naot, Shavit, and Wolfshtein 
(ref. 5). Their parameters, a , 3 , and y are related to 

these by 

a = GUA + GUB 

3 - GUB 

Y = 3 • GUC 

(It should be pointed out that the relationship between a , 3 , 

and y which is given in reference 5 is rescinded in reference 6 . ) 

One reason for not taking literally the process leading to equations 
(26) and (27) is that equation (23) does not agree with (18) when 
j and m are contracted. The fact that the latter is used in a 
term multiplied by the viscosity serves as one justification for 
the discrepancy. The use of equation (23) only as a suggestion for 
arriving at (27) is another. 


The substitution of the model given by equations (14) , (18), (26), 



< u i u pt + 5 


i j,£ 






- VUU 




- - pro [ A 2 (r 5 V u i“ ,n \j 1 + 
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3.- THE SCALE TENSOR 


Understanding of is sought through the scale tensor, and the 
scale tensor is defined in terms of the two-point velocity 
correlation. Consideration of two points a finite distance apart 
is cumbersome in general coordinates. Therefore, Cartesian 
coordinates are employed for most of this section. 

Consider two points in the same flow: - (x x , x 2 , x 3 ) and 

5 ~ (C. ■ » 5 » • C ) * Proceeding much as in the development of 
equation (12) , rewrite equation (11) with j substituted for i , 
interpret it as applying at ( , and multiply by ui (x) . Add the 

result to the original equation (11) interpreted as applying at x 
and multiplied by uj(£). Note that a function of x behaves as 
a constant in differentiation with respect to £ and vice versa. 
Using equations (8) and (10) and averaging, the result is 


3u i ( x)u j ( ° + . s * {£) au i< x)u j(0 + ^ 


It 




3x A 


. 0 tr. » c ( 


7 l 3u.(C) 


+ u!(x)u ,j "(c) — 1_IL + u ' £ ( x ) u ; (£) _i 
A 3?* J 


3u. (x) 


H 


3x 


.1 


3u’*(C)u! (x)u' (5) 3u ,A (x)uV. (x)u! (£) 

z J a- J 


H 


3x 


3P'(Ou[(x) 3P'(x)u'(0 


3£ J 


3x 


,2-r 


^ (O 3 2 u! (x)u! (£) V > - t'-o* 

+ v — „ — «• + r._. . . ... J 


o c 

H dV 


3x £ 3x £ 


(29) 


It is not obvious that this equation is equivalent to equation 
(12) in the limit as x and 5 approach the same point. The 
demonstration can be accomplished using techniques that are 
employed here for a different purpose. 

Sandri introduced the concept of modeling at this stage (ref. 7). 
One advantage is that the modeled equation should reduce to equa- 
tion (28) as x and y approach each other. This helps estab- 
lish the form of the models and also determines some of the 
constants. Thus, part of the arbitrariness in other approaches 
(see, for example, ref. 8) is avoided. 
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Another advantage is that it is possible to relate the properties 
of the modeled equation to the spectral properties of the turbu- 
lence of which so much is known. This aspect is not, however, 
pursued in this paper except to note that it has influenced the 
models. A third advantage is that the terms to be modeled have 
known properties that can be exploited in the construction of the 
model. 

Some technical apparatus is needed for what follows. New indepen- 
dent variables, c ** (t 1 ., t 2 , c 3 ), the centroid, and 
r * (r 1 , r 2 , r s ), the relative separation, are defined by 


Note that 


C = \ (x + X) 


r ’ » i x 


.13 

’ 7 ar 1 


a_ + 1 


a? 1 sx 1 a$ x 


13 1 3 

2 — x + 1 — i 
4 3X 1 4 3£ x 


A consequence of equation (34) is that for a function of x only, 
say 0 


, ' » - T , .. e * e(x) 

3x 3c x 

and for a function of £ only, say x 

. X - XCO 

3 r 3c 


it// : 




N ij ” ~ij 


0 ^<4 

R_. , + r* — 4 1 


3r 


The terms of equation (29) are now considered in groups. The first 

a. n Ot 


three, in terms of R. . , are 

_il + u & ( 5 ) — 11 + u £ (x) —|1 
35 3x £ 


3R. 


3 t^* + | (u*(C) + u £ (x>) — |l + (G £ (0 - u £ (x)) -4J 


3R 3 . , 

3r* 


(43) 

by equations (32) and (33). On the right-hand side of (43) and in 
similar circumstances below, £ and r are the independent 
variables; x and £ are merely abbreviations for £ ~ %r and 
C + %r , respectively. The next two terms, the "production terms," 
are • 


R i 


U l 


£ 3u.(x) » 8u-(C) j 3u-(x) 

. i « R * L_ 4 . R* ± 

3 3x £ 1 3£ £ j 3£ £ 


(44) 


by equations (36) and (37). 

The triple correlation terms are modeled: 


3 u ,£ (Du!(x)u!U) 3 u £ (x)u! (x)u! (£) 

— 5 i 4 + i 1 

3 r<: ■« ° ; -V c - «>>;■/• 


dX 


% 




;t\: fJ ° o ? . ’/o t tf. o 


« - vuu 


..85 


/ 3R. . 3R. £ 3R £ . 

qA ( _^l + _4_ + -a 
V 3C 3C J 3C , 


j 2vATC x 

+ - — x — (R. . - N..) - 


+ 2 BTC |[s (E.. -N ij)+ 4l Si . ( R \ - S\j 


(45) 
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(Here and in what follows, q and A and other functions of a 
single vector variable are to be evaluated at £ if not otherwise 
indicated.) Thus, it is hypothesized that these terms have a 
diffusive effect and also a pseudo-dissipative effect — "pseudo" 
since the last terms tend to drive Rjj toward , or away, 

depending on the signs of the constants, instead of toward 0. 
Turning to the terms involving pressure fluctuations, note that 


3P'(C)u[(x) 


by equation (37) . But 


3P'(5> 

3? 


u[(x) 


3P’U) 


H 


j 


u[(x) 


(46) 


3P’(5)u^(x) 3P'(5) 


3u! (x) 




— *— u[(x) + P*(0 — f- 
3C J 


(47) 


so, using equation (47) with (46) and making similar manipulations 
with the other P‘ term. 


3P* (v)u! (x) DP'(x)u'(0 


i 


3£ J 


ax' 1 


3P'(C)u!(x) 3P’(x)u:(0 


a 3 

The first pair is modeled: 0 ' :U 


1_ ru) 


3u^(x) 

3 ?" 


- p*(x) — i 


3u!(5) 


3C 
(48) 


3P’(5)u!(x) 3P‘(x)u!(0 

, i + , J 


35- 


35 


PMU 


+ PMU 


35 J L 
a 


A 2 l — i- + 


3R.. n l au i (x) aR*’”' "• 


1 3 ^n 


m. ' 

cZ 


4 a c n 


35 


35 


A 2/l 3 u *(Q + 1 S ^ (g) ^ 

a; 11 a ^ 4 35 n 3?^i 
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•1/ * 


+ PUW 


[ «3. ( 


3R*\ . / ™ . 

qA — r- 1 + — r I jjf 

3C / K X 3 « 


•onr 9. 3A 3 A /« . n ) (49) 

- ?0C t~ J ^ (R ij iy 

Here, in addition to the two types of diffusion terms foundin 
equation (28) , another term is included. Agaxn. this term unds to 
drive R.l toward (or away from) N., but o nly where t here xs a 

gradient^ in A . It is known (ref. 9) that P'U)u[(3c) is zero xn 

isotropic turbulence. Therefore, all the terns on right-hand 
Xde il equXion (49) must be zero iothetsituatron. in the last 
term that is assured by the factor (3A/3; ) • (3A/K ) • 

The other pair of terms on the right-hand side of equation (48) xs 
modeled: 

3u! (x) 3 uT(C) 

p-«) -4— +■»■<*) ri 

3C J 

r ^ «■) „ 5 "i« n 

- GUA hi— R i its R 3 


1 (^S£i r n +2j£&sl r" 

' 3 g ijU 5 n 1 3t n 


* Gb2 | [ h} - S ±j - | Sy (R\ - *“„)] < 5 °> 

^ItV 0 ? Sde Sat 6 is the 

Sfs^for’tfe IntSSnV S tj in the last term (according to 
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Property II). The other terms satisfy the requirement by virtue of 
t he fact that homogeneous turbulence will not exist in the presence 
of gradients of the mean velocity; 

For the final terms in equation (29), observe that 


3 Z u*(x)u!(C) 3 2 u!(x) 

— JL— — as — - — Ui 


ui(C) 


35 35 

or, using equations (36) and (37) 


35*35* 


3u! (x) 

> J- 

3u’(C) 

* 35* 

35* 


3 2 u!(C) 


+ U (x) 


" i ''" 35 V 


3x*3x* 


35 


35 


3 2 u^(s 

3S*3£* 


Thus, the last two terms of equation (29) can be written 

.2, 


/ 3 2 u! (x)u! (C) 3 2 u[(x)u' U)\ 3 2 R.. 

v — + o o — J — - v —rr 


3£*3£* 


3x*3x* 


35*35* 


2v 


3u! (x) 3u!(£) 


35 J 


35 


(51) 


One more model is needed: 


3ul(x) 3u- (g) 


35 


35 


g + ks. 
71 R ij + Av 


_a 

A' 


(• E ij + T- 4 8 ij E **) 


(52) 


c. -®S; ? y° ; c - c- ; >. V V* r- * " ,v - ■ **> ' >» c ' ' * 

Putting together the results of equations (43) through (52) , 
equation (29) becomes the modeled two-point correlation equation: 


3R. , 
3t 


-ii + 2-(u*(?) + u £ (x)) + - ^ £ < x )) ^r 1 


+ R. 


y 3u.(£) „ 3 u 4 (x) 

* J + R* . 


j 


x ' - vuu 


35" J 35 

+ 2v^C (R_ -N^) +2BTC | 


r' /«,, 3R * 9R*> 

— r i*( — r L+ — r + “T L 

35 b \35 35 J 3C J 


) - 
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* ~ f )9 
it if* 


- PMU 


SrT 


' 2/1 3u & < 

[ K [*7f 


3R. 


^(x) »"1 . 1 8n i (i ° SR“\1 

3C J L Si - * Sc 1 1 3c n 3C* jj 







^ 3£ n 

r. Tt & 1 


8C 


1 


3A 3A 



V 


T 3U|(f) 3Ui<x) 

+ GUA| — — — R- +— 5 R i 

L H n 1 3C n 3 

- U 3 h n + *\)] 

+ GCBf (*£L - h * 4^ - *,1 

_\3C C H 2 ) \3C H 1 J J 

a n" /• _\ % r 


+ GUC 


V 3S 1 


+ v 


35 

,5? - 2 ? R u - 2b 1 K + 


<0 + — R" n + GU2 fL, - S y - 1 q 13 (R” - S n J 


3? J 


/ 


3-3 

1 


g. . R 
& ij 


-) 


(53) 


It is easy to show that equation (53) reduces to equation (28) as 
x c goes to zero. The terms containing Njj drop out by equation 

(41) and S-^ drops out by equation (39). The term ^u (£) 

- u^(x^3R£j/3r^^ becomes zero since £ and x both approach £. 

The rest reduces term by term to equation (28) . 


The scale tensor, A^j , is defined through the quantity 


n ij 



hjrF 


dr dr dr~ 


(54) 
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(It is seen that 8^ is symmetrical in i and j .) Then, 


A iJ 


hi 

q 


(55) 


Actually, most of the analysis is done in terms of • Accord- 
ingly, equation (53) is divided through by 4r | r i 2 mi integrated 
over all r space. Many of the terms are expressible immediately 


in terms of 


a ij 


For example. 


i 



1 3 / 

PUW Z-j qA —j~ 

h 3 \ sr 


,1,2,3 
dr dr dr 

Mr? 


PUW 


= PUW 


3C 

3 

3C 





R i £ dr 1 dr 2 dr 3 
4w | r | 2 


since C and r are independent and^ q and A are, 
convention adopted above, functions of C only. Itis -een by 
equation (40) that S ii drops out and by equation (42) that 
drops out. This is in 3 contrast with the limit as r goes to 
zero where the whole terms containing dropped out. 

Terras containing factors evaluated at x or ? need an additional 
assumption before the integration can be completed. Roughly 
speaking, the assiimption is that the variation of Rij wtk C 
is smoother than the variation with r , that is, the flow^i^ 
locally homogeneous. More formally , the assumption is that it is 
sufficient to use the leading term in a moment expansion. 
Explicitly, it is assumed that y r ■■■* ■ *' * .* / ' 


R ± j (c,r) * ( t) 4 tt j r | 6 (r) 


(56) 


where 5(r) is the three-dimensional Dirac delta function defined 
by 



<j.(r) <5(r) dr 1 dr 2 dr 3 » 4>(0) 


(57) 


for any continuous function <|> 
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*7, 


?i 







or* by equations (8) and (35) 


Q. . - 0 


The equation for ft., can now be written. Since the r depen- 




(59) 


dence has been eliminated by the integration, only one point is 
involved in the equation so general coordinate notation is employed 
once more. The contraction (trace) of ft^ is written 

a - a\ (58) 

To take advantage of combinations of terms that occur, write 

AOM = a + ATC 
BOH = b + BTC 

The equation is 

“i) t + °ij.* + + “‘j 

. VDlj[qA(a 1J - t + 0* tt j +8 t j>t )].t - 

" - ^“(^(z 5 \n + i 5 i.n 

- P>»[a 2 (| u * >0 + l 5 j;n """.a)].! ' 

- PUW [(qA + (qA 8^';,) >t ] 


+ POC | A-*A >t Sly 


+ GUk ( Uj n + u in !)”j - | gjj u E in Sl n j) 
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f ' 



+ GUB [(%.n ' VjKl + (“i.n - %.i) a ”j] 
+ guc (u^ + S^) o + GU2 J(a tj - § gyfl) 


+ Vft. . 


A 


2 


vAOM 




2BOM (60) 


From this equation an equation for A.. can be derived, as 
follows. The contraction of equation J (28) (displayed in the next 
section) is an equation for q 2 . If it is multiplied by A- • and 

J 

subtracted from equation (60), with q 2 A. . written for a.. , the 

1J 

result, after dividing through by q 2 , is an equation for Ay . 

It is considerably longer than (60) and not worth writing out. To 

determine Aj . , it is simpler to solve for fi. . and divide the 
I* J V J 

result by q z . 

Thus, it is equation (60) that has been programmed for solution, 
along with equations (8), (9), and (28). Actually, for the model- 
ing parameters that appear in both ( 28 ) and ( 60 ), the program 
version of ( 60 ) uses different names. In addition, the program 
version has two extra terms to allow direct comparison with earlier 
versions of the scale equation. Details are given in reference 1. 
In all the runs made so far, the renamed parameters have been 
given the same values as their prototypes, and the coefficients of 
the extra terms have been set to zero. 


4.- THE BOUNDARY LAYER ASSUMPTIONS 


In this initial attempt to study the behavior of the system 
consisting of equations (8), (9), (28), and (60), two-dimensional 
steady flow over a smooth fiat plate is chosen. A Cartesian coordi 
nate system, with C 2 , C 3 denoted by x, y, z is used, taking 

the free stream parallel to the plate in the x direction and the 
y axis perpendicular to the plate so that the flow is envisioned 
in the x-y plane. Derivatives with respect to z , as well as to 
time t , are zero. The velocity components u lt u 2 . u 3 _ are 
denoted u, v, w . By the two-dimensional assumption, w = 0 , 
but this doesn’t apply to w' . 

Using x and y as subscripts to denote partial derivatives with 
respect to x and y , equation (8) becomes 

^ + v y = 0 (61) 

Equation (9) , for i = 1 , becomes 


uu x + vu y + (u'u') x + (u* v* ) y - - P x + v(u xx + u yy ) (62) 

and for i -2 , 


uv x + vv y + (u'v’) x + (v’v') y - - Py + vCv^ + Vyy) (63) 

Since the flow is two-dimensional, i = 3 in equation (9) is not 
relevant. At this point, the boundary layer assumptions are 
invoked. 

■'••• . > C ^ w o V V ' ° •' ° . :••••. h.. ■ c ..**•' o y <-• ' ... • * 

At a high enough Reynolds number and in the absence of extremes in 
the pressure gradient, it is generally recognized that the 
influence of a flat plate on the flow along it is confined to a 
narrow boundary layer with the following characteristics. Deriva- 
tives with respect to y are large compared with derivatives with 
respect to x , v is small compared to u , and the two effects 
balance each other so that the two terms in the continuity equation 
(61) are equal in magnitude. 

In equation (62) then, u X x is negligible compared to uyy and, 

if it is assumed, that the various components of are of the same 

order of magnitude, (u'u') x is negligible compared to (u'v')y 
On the other hand, the first two terms are comparable. Therefore, 
equation (62) can be approximated by 
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uu + -vu 
x y 


' P x - 


(u’v')„ + V 


Uyy 


(64) 


where (u’v')y has been put on the other side of the equation 
since it is interprets d_as representing a stress contributing to 
the rate of change of u along a streamline; that is, the left side. 

The process of arriving at equation (64) from (62) can be formalized. 
Interpret the equations as being in terms of nondimens ional variables 
so that v stands for the reciprocal of the Reynolds number. Let 
6 be the nondimensional measure of a nominal boundary layer thick- 
ness, assumed small. Then u and v are taken to be of order 6° 
and 6 1 , respectively; the operation of differentiation with 
respect to x and y is taken to be of order 6° and 6" 1 , 
respectively; v is taken to be of order S 2 . All this is standard 
in the theory of laminar boundary layers. The components of 

u!u! are taken to be of order 5 l . The terms that were dropped 

from equation (62) ere, by these criteria, of order 6 l or 6 2 
and hence small compared to those found to be of order 6° ; that 

is, small compared to the rest of the terms except, perhaps, P x 
which is considered in the next paragraph. 


or 


The terms of equation (63) are seen to be all of order i 

smaller except (v'v')y which is of order 6° and Py which must 
also be of order <5® 


Thus. 


<vv) y - - r y 


which can be integrated to give 


v* v 1 


(vV), 


- P + P 


where the subscript stands for free stream conditions (i.e., the 
edge of the boundary layer) . Solving for P and differentiating 
with respect to x gives 


P v - P e - (W) + (vV), 

x e x x e x 


for use in equation (64). However, the last two terms are of order 
6 1 and so are dropped. From now on, F x will be taken to mean 
, the external pressure gradient. 


The limitation on the pressure gradient needed to justify the 
boundary layer assumptions can now be made less vague by stating 
that P x must be at most of order 6° .It may, of course, be 
smaller, in which case it won't contribute appreciably to equation 
(64). 
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Similarly, it should not be assumed that the other terms in (64) are 
always of the order assigned to them. Indeed it is well known that 
in the main part of a turbulent boundary layer the term vuyy is 

negligible compared -(u'v')y but that the reverse is true in the 
sublayer next to the wall. 

The sane technique is applied to equation (28) . The modeling 
constants are assumed to be of order 6° . The definition of q 2 
can be rewritten 


2 — 

q *» u'u' + v'v' + w'w’ 

. h 

so q is of order S' 5 . To estimate the order of A , consider 
the contraction of equation (28): 


(q 2 ) t + u £ (q 2 > f * + 2u’V n ^ - vuu qd((q 2 )‘* + 2(u' V n ) ^ , 




2WJw[qMu , V n ) jJl ] jn + v(qV A t 


? Va n 2 

A r 


2b | q 2 


(65) 


Among the urnno deled terms, the dominant one, after mak ing the 
boundary layer assumptions, is a "production term," 2u r v r Uy , of 
order 5° . It is generally recognized_that this term (which is 
ordinarily negative, since u T v T " and Uy are usually of opposite 
sign) is principally balanced by "dissipation," that is, the last 
two terms in equation (65) . - If A is assigned the order 6' with 

Y to be determined, the first of the dissipation terms 2(va/A 2 )q 2 

has order 6 3-2 Y and the second 2b(q/A)q 2 has order 6 3/2 ~ Y 
For either of these to balance the production term of order 6° , 

Y * 3/2 . The conclusion then is that taking A to be of order 
6 3,2 gives the proper ordering. 

The order of all the terms in equation (65) can now be determined. 
Keeping only the largest gives 


2u‘ v* u = - 2 q 2 - 2b f- q 2 
y ^ a 


( 66 ) 


This is a useful approximation for many purposes and is exploited 


it' * 
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in the next section. (As is the case with (u'v')y and vuyy in 
equation (64) , the two terms on the right of (66) are not important 
in the same part of the flow, Indeed, their ratio is 


k ak „ 

a v 



which varies from zero at the wall to values large compared to one.) 


However, equation (66) is too crude an approximation for a study of 
the evolution of the boundary layer, so terms of order 6 1 are 
retained. Equation (65) becomes 


u(q 2 ) x + v(q 2 ) y + 2u'u‘ u x + 2u’v' u y + 2v'v' v y 
- VUu[qA((q 2 ) y + 2(VV) y )] y 

» - 2PUW^qA(v'v') y J y + v(q 2 )yy " 2 ^ ~ 2b ^ ^ (67) 


The changes in the modeling that have been introduced since refer- 
ence 4 all drop out in the contraction of equation (28). Therefore, 
equation (67) is essentially the same as we would have written it 
then, with the exception of the terms 2u*u' u x and 2v T v r Vy 
which we would have dropped, as do other workers in the field. The 
justification, if anything is said about it at all, it usually 
something like: "They are production_terms ; they are small compared 

to the other production term (2if v’ r Uy) ; therefore they are 
negligible.” But if they are negligible, so are other terms in 
equation (67) ; the result is (66) again. 

Note that insofar as these considerations apply to a comparison of 
the "small" production terms (which are also known as the stream- ° 
tube stretching terms) with the convection terms, u(q 2 ) x and 
v(q 2 ) v , the reasoning is independent of the order of magnitude 
assigned to A , or for that matter to u^uj . In fact, it is hard 

to conceive of a systematic ordering that would drop those produc- 
tion terms but retain the convection terms. 

The boundary layer approximations can now be appl ied to equa tion 
(28) without contraction. It is found that u w and v w form 

an independent set and do not appear in the other equations; there- 
fore they are ignored. So far, the equations have been written with 
the production terms and the triple correlation terms on the left 
with the convection terms, to emphasize their common origin in the 

term (u*u.) of equation (3). Henceforth, as is customary, only 

X * Xf 


28 


I5<' 



i 


the convection terms are written on the left, in parallel with 
equation (64). 

In the boundary layer approximation, the terms multiplied by PUW 
appear in the same combinations as certain of the terms multiplied 
by VUU. Therefore it is convenient to introduce a new parameter: 

VPW - VUU - PUW 

The final version of equation (28), then, consists of four equations: 
u(^) x + v(uV) y = -2^V u x - 2^ u y + VUu[qA(^ T u r ) y ] y 

+ | GUA.(2u'u' u x + 2u'v' u y - v'v' v y )+2 GUB u'v' u y 

2 

+ 2 GUC q 2 u x + GU2 |(u’u’ - ^) + v(u ’ u,) yy 
- 2 uV 2b ^8u'u' + 1 y — q 2 ) (68) 

u(^v r ) x + v(v r v r ) y - -2VV v y + VUU qA(v r v*) y jy + 

+ 2 VPW qA(v r v r )y]y 

* | GUA(uV u x + uV u y - 2vV v y > -2 GUB u'v' u y 

2 ] 

+ 2 GUC q 2 v y + GU2 |(v'v T - ^) + v^'v')^ 

_ 2 ^ y? - 2b J(&v’^ + q 2 ) (69) 
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u(w'w') x + v(w e w') y - VUu[qA(w*’w , ) y ] y 

- | GUACu'u* u x + uV Uy + v'v' Vy) +G02 ^w'w' - 

+ vCw’w'Jyy - 2 ^ wV - 2b J^w'w’ + q 2 ) (70) 

u(u'v') x + v(u‘v')y ■ - v'v 1 u y + VUU^qACu'v’JyJy 
+ :VFw[qA(tt^) y ] y - | PMu[A 2 B y (v'v’) y ]y 
+ GUA v'v' Uy - GUB(u'u* - v'v‘)u y 

+ GOO q 2 Uy + GU2 | u'v' + v(uV)yy 

- 2 ~| uV - 2b jsu’v' ■' (71) 

A 2 A 

Applying the same approximations to equation (60) for. fly involves 
no new concepts or problems. The smallest terms retained are again 
of the order of the convection terms and, hence, of Gfj itself, 

namely, 6 5/2 . As is the case for u'w' and v’w* , fii 3 and 
Slzz do not impinge on the other components and are ignored. 

Some new notation is introduced, largely for convenience in 
programming: 

ee ** 

FF " Q 22 
EF - fl 12 

where each letter pair is taken to be a single symbol; also 

VOW « VOM - POW 
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The final version of equation (60) „ then, consists of these 
equations : 

5EE x + vEE y » -2EE u x - 2EF Uy + VUU(qAEE y ) y 
+ ®0C J A y 2 EE 

+ | GUA(2EE u x + 2EF G y - FF v y ) + 2GUB EF u y 
+ 2GUC a a x + GU2 J. ( EE - |) + vE Eyy 

- 2 EE - 2 BOM EE + (72) 

5FF X + vFF y * -2FF v + VUU(qAFF y ) v + 2VPW(qAFF ) 

.■ ^ y y y y 

+ POC 3. A y 2 FF - | GUACEE .5^.+ EF u y - 2FF v y ) 

- 2GUB EF u y + 2 GUC Si v y + GU2 | ^ FF - 

■ + v FF yy 2 FF - 2 BOM |{VfF -r fl) (73) 

uGG x + vGG y = VUU(qAGG y ) y + POC J A y 2 GG 

- | GUA(EE u x + EF u y + FF V y ) + GU2 J ^GG - 

+ v GG yy - 2 GG - 2 BOM 9-^BGG + (74) 
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uEF x + vEFy “ - FF u y + VUU(qAFF y ) y + VFW(qAFF y ) y 


- \ PMU(A 2 u y FF y ) y + PQC J A y 2 EF + GUA FF G y 

- GUB (EE - FF)Uy + GUG Q u y + GU2 J EF 

+ V EFyy - 2 EF - 2BOM J 0 EF , (75) 





5.- RESULTS 


5.1 The Program 

r6A? n?> ^ r /’Ifl\ 8 f^ ia t0 u S /^? the set consisting of equations (61), 
(6A) and (68) through (75) by finite-difference methods is 
described in reference 1. A preliminary series of runs has been 

?60? iSnSXlfr 8 *? 16 S^ tU j e of . th< r ®cale tensor insofar as equation 
i?n2*.SF* It xt l 1110 description of the program in this paper is 
limited to those features actually used in the runs discussed. 

o f the runs had a number of features in common. These include 
boundary conditions: all the dependent variables are 

f - he ^f 1 }’ . y • 0 * and all go to zero in the free stream 
\ which ?- s set to a constant, u- , there and v which 
takes whatever values are needed to satisfy continuity The 
pressure gradient P x is set to zero. 

A for use in the models at 


The program computes two versions of 
each point, A r and Ag . (The subscripts are m 
ratio and empirical.") The first is defined by 


that is, 
function 


* n 

A r is the contraction of A 


ij 


The other is the 


A„ = min 
e 


c 6 


99 


d y 


T^ere _c and d are constants and 6 gq is the value of 
which u .99 Ug . The values used for the constants are 


y for 


c — .17 

d - .65 

which, in conjunction with the modeling parameters given below, are 
known to give good results for a flat plate boundary layer when A 
is used for A . Either A r or A e was used throughout any 8 
one run. e 

5.2 Evaluation of the Paramet ers 

^deli ng parameters appearing in the boundary layer equations 

^ or u i u j ' e< ? ua £ions (68) through (71), were given the following 
values : 
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VUU «■ .3 
VPW - 0 
PMU “0 
GUA ■ 0 
GUB - 0 
GUC ^ 0 
GU2 » -1. 
a “ 3.25 
b « .125 
6*0 

(Hence, the extra parameters of reference 1, referred to at the end 
of section 3 , were given the values : VOM = . 3 , VOW = POM « GOA = 

GOB ** GOC * 0 , G02 * -1 . , SAM ■ TAL =0 . ) These parameters gave 
good results when used with A e for A for the mean flow and for 

u£ui as they had in earlier programs. One exception to that state- 
ment is that the calculated maximum values of u'u' , v’v* , and w'w' 
are found to be approximately in the ratios 


u'u* :v'v' :w’w’ - 2:1:1 

which is not in agreement with experiment. This defect can easily 
be remedied by changing the values of GUA, GUB, and GUC, but that 
was not deemed necessary for this preliminary study, especially 
since various experiments do not agree on what the ratios should be. 

There remains a group of three parameters, AOM, BOM, and POC, for 
which values are needed. Consider equation (67) in the free stream 
where u is constant and derivatives with respect to y are 
negligible: 

u(q^) x “ -2 ^ - 2b (76) 

The corresponding equation for 0 is obtained from the sum of 
equations (72), (73), and (74) (since 0 * EE + FF + GG) with the 
same assumptions. It is 

uft = - -4° -- fi - 2 BOM S- 0 (77) 

X A 
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For this discussion, A is identified with A r , so substituting 
q 2 A for ft gives 


u A <q 2 ) x + uq 2 A x « - q 2 A - 230H J q 2 A 

Subtract equation (76) multiplied through by A from this and then 
divide the result through by q 2 : 


u A x » - a - 2 (BOM - b) J A = - — ^ - 23TC q (78) 

A 

For high Reynolds numbers. Re. « qA/v , the terms involving v 
can be neglected. Then, 


■ i> 


i(q 2 ) x - -2b V 


u A x « -2BTC q 


(79) 


(80) 


It is known (ref. 10) that in these circumstances q 2 decays like 

r n where n is probably around. 1.25. For this to be compatible 
with equations (79) and (80), the. relation 


BTC 


(z ' n) 


must hold. Thus , 


BOM - b + BTC 


(z ~ n) 


b 


For n ® 1.25 and b * .125 , these relations give BTC = -.0375 
and BOM = .0875 . As expected, A x is positive so that A grows 
•as. q 2 decays. 

While it is possible to approach the evaluation of AOM by 
considering equations (76) and (77) in the limit of low Reynolds 
numbers, a different flow situation is used here. It is easy to 
see that in the limit as y goes to zero, the dominant terms of 
equation (67) are v(q 2 )yy - 2(va/A 2 )q z so 


<q 2 > 


yy 


- 2a %r 


( 81 ) 


35 


■:y/f 


very near the wall. The corresponding equation for fl is 

ftyy = 2A0M (82) 

With A (still identified with A r ) given by ft/ q 2 , this pair of 
equations can be solved by assuming that q 2 varies like y? and 
ft varies like y H . It is easily established that q = p + 1 ; 
that is, A is linear with y , a result in good agreement with 
views of the subject since Prandtl. If the slope of A as a 
function of y , for small y , is designated d (as it is for A,,) 
it is found further that p(p - 1) = (2a/d 2 ) or, taking the posi- ’ 
tive root, 

8a"\ 

?/ 

Also, 

AOH _ q(q. - 1) _ p + 1 
a pCp -T) p -~ T 

For a > 3.25 and d = .65 , p = 4.454 and AOM - 5.132 , 
approximately. 

To obtain a value for POC , consider equations (68) through (71) 
keeping only the highesfc-order terms (as in eq. (76)) with the 
additional assumption that Re^ is large. The resulting set is 
equivalent, insofar as the modeling is the same, to the "superequi- 
librium" equations of reference 11. It is 


(2 - | GUA - 2GUB) u'v' u * 

GU2 ^ ^u'u' - - 2b ^8u'u' + — q 2 ^ (83) 


(| GUA + 2GUB) u*v‘ u y = 

■ Gu2 J “ %-j - 2b l^v'y 1 + q 2 ) (84) 
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| GUh uV u y - GU 2 }(wV - 4 ) - 2b K BW ' W ‘ +? ^ q2 ) (85> 


- GUA - GUBj v*V u y 


-GUB u'u’ u y + GUC q u y 


2 u v + GU2 ^ u'v’ - 2b ^8 u'v’ 

( 86 ) 


These are regarded as 

K!n| £ arVdue%i ?SeS ? ( $ ba ; a “ “ ™ dlly 

obtained. Add equations (83), (84), and t ) 


2u' v' Uy * *2b ^ q 


(87) 


,,,, v q«bstitutinR this result back inM equations 
: c °tnpare eq. (66).) Substitut g so lved for u -rjjr , v v 

It S S| y th°at Se eonSls 

squations (84) s|d_ J.SW^.t «y st aJlpMrin ,ttoa w If these 
whether or not v v f rf . - n J; , «.*,* value of i?v^ from 

iS? (M>. aha — u can be " rl “ en 

MOW - 2 b 6 ) 2 q 2 - | Cl - SDA) [GU2 + »» * »>1 

+ 2bf(l - GUft)(i G8A +2 GUb) - 2 GBB 2 ] 


(88) 


+ (GU2 - 2b8)GUcl A 2 u y 2 

or the values of the parameters given above, this becomes simply 

- 2 (89) 


= 2 


aV. 


£lie^ dJtt^ 

.ations . 
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It is laaown (refs. 8 and II) that this approximation represents the 
situation in the law - o f - th e - wa 1 1 region of a turbulent boundary 
layer. For example, equation (88) is in agreement with the known 
properties of that region in that the level of turbulence is 
roughly constant, A is proportional to y .and the velocity 
profile is logarithmic so that u^ is proportional to y~ l , 

A similar analysis can be carried out on equations (72) through (75) 
but, in view of the variation of A in the law-of-the-wall region, 
terms involving A or are retained, and it is assumed that 

/ J y 

each component of ^ is proportional to A^ . For the values 

of the parameters already mentioned, the equation corresponding to 
equation (89), to skip a few details, is 


A*Uy 2 [.175 - (.3 + p °C)A y 2 ][l - (.3 + P0C)A y 2 ] 2 

or, using equation (89) 

(.175 - z)(l - z) 2 = .275 

where z ~ (.3 + ?QC) A^ 2 . The real root of this cubic is approxi- 

mately z - -.066&9 . If Ay is assumed to have the value it has 
for A - A e , namely, A y - -.65 , it is found that POC - -.4578 . 

To summarize, the following values have been found for the new 
parameters: 

AOM - 5.132 
BOM = .0875 
POC * -.4578 

These, of course, are to be considered just first guesses with 
adjustments to be made after trial runs. 

5 . 3 The Runs 

A preliminary run was made to establish initial conditions for the 
other runs. The initial conditions for the preliminary rum were 
those of a laminar boundary layer with a spot of low-level turbu- 
lence. Hie initial conditions for flj* and the values cf AOM, 
BOM, and POC were arbitrary. The rum used A = A e ; it was stopped 
at a Reynolds number based on x , Re v = (u x/v) , of 5 million/ 

X w 


it jr. 
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After verifying that the profiles of u » v , and u|uj were 

representative of a fully developed turbulent boundary layer at that 
Reynolds number, they were used as initial conditions for the rest 
of the runs. Initial values of ftij were set according to the 
formula 



A e°I u i u j + A e (1 “ 


a 


I 



where Or was either ,9 or 1 . For any value of , this 
formula guarantees that the A r computed for these initial condi- 
tions is equal to . 

The first such run used the values of the new constants, ACM , BOM , 
and POC , found above, had oj = 1 and used A * A r . As the 
solution proceeded downstream, A became quite small through most 
of the boundary layer. This emphasized the dissipation terms so 
that the level of turbulence, and with it the skin friction, dropped 
off very rapidly. A few trials with modified values of the new 
parameters showed similar behavior ox the opposite, with. the scale 
and turbulence levels growing much too large. It was quickly 
apparent that it would not be easy to find suitable values. 

A series of rimts was made with A = so that the turbulence level 
would remain normal while the effect of a variation of the new 
parameters on fci and on A r could be monitored* A set or values 
that gave good results in this mode was 

AOM ~ 5.65 

BOM - .0875 
POC = -.27 

However, with A *» A r , the solution again showed rapidly increasing 
A and turbulence levels. Systematic changes in the parameters 
were tried to no avail. The results were unstable in the sense that 
a small change in one of the parameters could produce a large change 
in the solution. This indicates that even if one satisfactory run 
were obtained, it would no doubt be sensitive to initial conditions 
and other factors that should be extraneous to a useful, calculation 
procedure. 

These results are discouraging but not disastrous. For one thing, 
only the new parameters were varied. It is not necessarily true 
that values of the other parameters good for A * A ? remain good -• 
for A « . For another thing, the new modeling is deficient in 

a way that^seems obvious in retrospect. A look at equation (53) 
sugpests that the "POC" term could be expanded in the pattern of the 
"ATC" and "BTC" terms. That is, the "POC" term could be replaced by 
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5.4 -The . Nature of the Tensor 

sSS S.S»'=J- , -S 


fl ij " I 


s ij n + s Lfe - 1 Sij q 2 ) 


(90) 


is sufficient to ToTeYj ’ lt 

Equation (90) was tested by forming the ratio 
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at each point for each pair i,j . If the various Sj^ at a 

given point are equal, the hypothesis is verified. It was found 
that the diagonal elements (i <* j) were very close to each other, 
generally within one part in a thousand, but the off-diagonal value 
tended to be about ten percent higher. 

A further hypothesis is that 



where a is a constant. This too was tested, using the value 
found on the diagonal of for Sl . It was found that o 

was almost always in the range from .85 to .95, and most often 
between .88 and .92. Except for some transients dependent on the 
initial conditions, these results concerning and a were 

remarkably constant across the rims, independent of whether A *= Aj- 
or A “ A~ was used and independent of the choice of parameters. 
Too much shouldn’t be made of these observations until the hypo- 
theses are tested in other flow situations. If they hold up, how- 
ever, it appears that if the accuracy of the first is acceptable, 
then setting a - .9 is equally valid and only one differential 
equation need be solved to determine A r , 
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